Analysis of blood and nasal epithelial transcriptomes to identify mechanisms associated with control of SARS-CoV-2 viral load in the upper respiratory tract (preprint)

Background: The amount of SARS-CoV-2 detected in the upper respiratory tract (URT viral load) is a key driver of transmission of infection. Current evidence suggests that mechanisms constraining URT viral load are different from those controlling lower respiratory tract viral load and disease severity. Understanding such mechanisms may help to develop treatments and vaccine strategies to reduce transmission. Combining mathematical modelling of URT viral load dynamics with transcriptome analyses we aimed to identify mechanisms controlling URT viral load. Methods: COVID-19 patients were recruited in Spain during the first wave of the pandemic. RNA sequencing of peripheral blood and targeted NanoString nCounter transcriptome analysis of nasal epithelium were performed and gene expression analysed in relation to paired URT viral load samples collected within 15 days of symptom onset. Proportions of major immune cells in blood were estimated from transcriptional data using computational differential estimation. Weighted correlation network analysis (adjusted for cell proportions) and fixed transcriptional repertoire analysis were used to identify associations with URT viral load, quantified as standard deviations (z-scores) from an expected trajectory over time. Results: Eighty-two subjects (50% female, median age 54 years (range 3-73)) with COVID-19 were recruited. Paired URT viral load samples were available for 16 blood transcriptome samples, and 17 respiratory epithelial transcriptome samples. Natural Killer (NK) cells were the only blood cell type significantly correlated with URT viral load z-scores (r = -0.62, P = 0.010). Twenty-four blood gene expression modules were significantly correlated with URT viral load z-score, the most significant being a module of genes connected around IFNA14 (Interferon Alpha-14) expression (r = -0.60, P = 1e-10). In fixed repertoire analysis, prostanoid-related gene expression was significantly associated with higher viral load. In nasal epithelium, only GNLY (granulysin) gene expression showed significant negative correlation with viral load. Conclusions: Correlations between the transcriptional host response and inter-individual variations in SARS-CoV-2 URT viral load, revealed many molecular mechanisms plausibly favouring or constraining viral load. Existing evidence corroborates many of these mechanisms, including likely roles for NK cells, granulysin, prostanoids and interferon alpha-14. Inhibition of prostanoid production, and administration of interferon alpha-14 may be attractive transmission-blocking interventions.

Safety and immunogenicity of a variant-adapted SARS-CoV-2 recombinant protein vaccine with AS03 adjuvant as a booster in adults primed with authorized vaccines (preprint)

Background Booster vaccines providing protection against emergent SARS-CoV-2 variants are needed. In an international phase 3 study, we evaluated booster vaccines containing prototype (D614) and/or Beta (B.1.351) variant recombinant spike proteins and AS03 adjuvant (CoV2 preS dTM-AS03). Methods Adults, primed 4-10 months earlier with mRNA (BNT162b2, mRNA-1273]), adenovirus-vectored (Ad26.CoV2.S, ChAdOx1nCoV-19) or adjuvanted protein (CoV2 preS dTM-AS03 [D614]) vaccines and stratified by age (18-55 and [≥]56 years), were boosted with monovalent (MV) D614 (5[≥]g, n=1285), MV (B.1351) (5g, n=707) or bivalent (BiV) (2.5[≥]g D614 plus 2.5[≥]g B.1.351, n=625) CoV2 preS dTM-AS03. SARS-CoV-2-naive adults (controls, n=479) received a primary series (two injections, 21 days apart) of CoV2 preS dTM-AS03 containing 10g D614. Antibodies to D614G, B.1.351 and Omicron BA.2 and BA.1 variants were evaluated using validated pseudovirus (lentivirus) neutralization (PsVN) assay. D614G or B.1.351 PsVN titers 14 days (D15) post-booster were compared with pre-booster (D1) titers in BNT162b2-primed participants (18-55 years old) and controls (D36), for each booster formulation (co-primary objectives). Safety was evaluated throughout the trial. Results of a planned interim analysis are presented. Results Among BNT162b2-primed adults (18-55 years old), PsVN titers against D614G or B.1.351 were significantly higher post-booster than anti-D614G titers post-primary vaccination in controls, for all booster formulations, with an anti-D614G GMT ratio (98.3% CI) of 2.16 (1.69; 2.75) for MV(D614), an anti-B.1.351 ratio of 1.96 (1.54; 2.50) for MV (B.1.351) and anti-D614G and anti-B.1.351 ratios of 2.34 (1.84; 2.96) and 1.39 (1.09; 1.77), respectively, for BiV. All booster formulations elicited cross-neutralizing antibodies against Omicron BA.2 across vaccine priming subgroups and against Omicron BA.1 (evaluated in BNT162b2-primed participants). Similar patterns in antibody responses were observed for participants aged [≥]56 years. No safety concerns were identified. Conclusion CoV2 preS dTM-AS03 boosters demonstrated acceptable safety and elicited robust neutralizing antibodies against multiple variants, regardless of priming vaccine. ClinicalTrials.gov: NCT04762680

Pseudotemporal whole blood transcriptional profiling of COVID-19 patients stratified by clinical severity reveals differences in immune responses and possible role of monoamine oxidase B (preprint)

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection is associated with highly variable clinical outcomes. Studying the temporal dynamics of host whole blood gene expression during SARS-CoV-2 infection can elucidate the biological processes that underlie these diverse clinical phenotypes. We employed a novel pseudotemporal approach using MaSigPro to model and compare the trajectories of whole blood transcriptomic responses in patients with mild, moderate and severe COVID-19 disease. We identified 5,267 genes significantly differentially expressed (SDE) over pseudotime and between severity groups and clustered these genes together based on pseudotemporal trends. Pathway analysis of these gene clusters revealed upregulation of multiple immune, coagulation, platelet and senescence pathways with increasing disease severity and downregulation of T cell, transcriptional and cellular metabolic pathways. The gene clusters exhibited differing pseudotemporal trends. Monoamine oxidase B was the top SDE gene, upregulated in severe>moderate>mild COVID-19 disease. This work provides new insights into the diversity of the host response to SARS-CoV-2 and disease severity and highlights the utility of pseudotemporal approaches in studying evolving immune responses to infectious diseases.

Unravelling the transcriptome of the human tuberculosis lesion and its clinical implications (preprint)

ABSTRACT Tuberculosis (TB) disease causes up to 1.5 million deaths every year and represents an important problem of public health at worldwide level. Here, we quantified the gene expression signatures of granuloma biopsies across human TB pulmonary lesions, and validated the best gene candidates using NanoString technology, profiling 157 samples from 40 TB patients who underwent surgery. We characterised the transcriptional profile of the TB granuloma in comparison to healthy tissue, described an 11-gene signature and measured 7 proteins in plasma associated with it. We demonstrated a gradient of immune-related transcript abundance across the granuloma substructure and evidenced metabolically-active Mycobacterium tuberculosis in the lesions. Patients who converted to sputum negative after two months of starting treatment, showed enriched inflammatory pathways in the lesion several months after, supporting use of sputum culture conversion (SCC) as a prognostic biomarker during clinical management and as a factor to prioritise patients when considering lung surgery.

A multi-tissue study of immune gene expression profiling highlights the key role of the nasal epithelium in COVID-19 severity (preprint)

Background: COVID-19 symptoms range from mild to severe illness; the cause for this differential response to infection remains unknown. Unravelling the immune mechanisms acting at different levels of the colonization process might be key to understand these differences. Methods and findings: We carried out a multi-tissue (nasal, buccal and blood; n = 156) gene expression analysis of immune-related genes from patients affected by different COVID-19 severities, and healthy controls through the nCounter technology. We then used a differential expression approach and pathways analysis to detect tissue specific immune severity signals in COVID-19 patients. Mild and asymptomatic cases showed a powerful innate antiviral response in nasal epithelium, characterized by activation of interferon (IFN) pathway and downstream cascades, successfully controlling the infection at local level. In contrast, weak macrophage/monocyte driven innate antiviral response and lack of IFN signalling activity were shown in severe cases. Consequently, oral mucosa from severe patients showed signals of viral activity, cell arresting and viral dissemination to the lower respiratory tract, which ultimately could explain the exacerbated innate immune response and impaired adaptative immune responses observed at systemic level. Results from saliva transcriptome suggest that the buccal cavity might play a key role in SARS-CoV-2 infection and dissemination in patients with worse prognosis. Conclusions: We found severity-related signatures in patient tissues mainly represented by genes involved in the innate immune system and cytokine/chemokine signalling. Local immune response could be key to determine the course of the systemic response and thus COVID-19 severity. Our findings provide a framework to investigate severity host gene biomarkers and pathways that might be relevant to diagnosis, prognosis, and therapy.

Characterisation of the blood RNA host response underpinning severity in COVID-19 patients (preprint)

Infection with SARS-CoV-2 has highly variable clinical manifestations, ranging from asymptomatic infection through to life-threatening disease. Host whole blood transcriptomics can offer unique insights into the biological processes underpinning infection and disease, as well as severity. We performed whole blood RNA Sequencing of individuals with varying degrees of COVID-19 severity. We used differential expression analysis and pathway enrichment analysis to explore how the blood transcriptome differs between individuals with mild, moderate, and severe COVID-19, performing pairwise comparisons between groups. Increasing COVID-19 severity was characterised by an abundance of inflammatory immune response genes and pathways, including many related to neutrophils and macrophages, in addition to an upregulation of immunoglobulin genes. Our insights into COVID-19 severity reveal the role of immune dysregulation in the progression to severe disease and highlight the need for further research exploring the interplay between SARS-CoV-2 and the inflammatory immune response.

Increased serum levels of sCD14 and sCD163 indicate a preponderant role for monocytes in COVID-19 immunopathology (preprint)

BackgroundEmerging evidence indicates a potential role for monocyte in COVID-19 immunopathology. We investigated two soluble markers of monocyte activation, sCD14 and sCD163, in covid19 patients with the aim of characterizing their potential role in monocyte-macrophage disease immunopathology. To the best of our knowledge, this is the first study of its kind. MethodsFifty-nine SARS-Cov-2 positive hospitalized patients, classified according to ICU or non-ICU admission requirement, were prospectively recruited and analyzed by ELISA for levels of sCD14 and sCD163, along with other laboratory parameters, and compared to a healthy control group. ResultssCD14 and sCD163 levels were significantly higher among COVID-19 patients, independently of ICU admission requirement, compared to the control group. We found a significant correlation between sCD14 levels and other inflammatory markers, particularly Interleukin-6, in the non-ICU patients group. sCD163 showed a moderate positive correlation with the time at sampling from admission, increasing its value over time, independently of severity group. ConclusionsMonocyte-macrophage activation markers are increased and correlate with other inflammatory markers in SARS-Cov-2 infection, in association to hospital admission. These data suggest a potentially preponderant role for monocyte-macrophage activation in the development of immunopathology of covid19 patients.

The impact of super-spreaders in COVID-19: mapping genome variation worldwide (preprint)

The human pathogen severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is responsible for the major pandemic of the 21st century. We analyzed >4,700 SARS-CoV-2 genomes and associated meta-data retrieved from public repositories. SARS-CoV-2 sequences have a high sequence identity (>99.9%), which drops to >96% when compared to bat coronavirus. We built a mutation-annotated reference SARS-CoV-2 phylogeny with two main macro-haplogroups, A and B, both of Asian origin, and >160 sub-branches representing virus strains of variable geographical origins worldwide, revealing a uniform mutation occurrence along branches that could complicate the design of future vaccines. The root of SARS-CoV-2 genomes locates at the Chinese haplogroup B1, with a TMRCA dating to 12 November 2019 - thus matching epidemiological records. Sub-haplogroup A2a originates in China and represents the major non-Asian outbreak. Multiple founder effect episodes, most likely associated with super-spreader hosts, explain COVID-19 pandemic to a large extent.